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ABSTRACT 

We present the discovery and follow-up observations of 142 pulsars found in the Parkes 
| 20-cm multibeam pulsar survey of the Galactic plane. These new discoveries bring the 

total number of pulsars found by the survey to 742. In addition to tabulating spin 
' and astrometric parameters, along with pulse width and flux density information, 

we present orbital characteristics for 13 binary pulsars which form part of the new 
q ■ sample. Combining these results from another recent Parkes multibeam survey at high 

Galactic latitudes, we have a sample of 1008 normal pulsars which we use to carry 
, out a determination of their Galactic distribution and birth rate. We infer a total 

d • Galactic population of 30000 ± 1100 potentially detectable pulsars (i.e. those beaming 

towards us) having 1.4-GHz luminosities above 0.1 mJy kpc 2 . Adopting the Tauris & 
i— i , Manchester beaming model, this translates to a total of 155000 ± 6000 active radio 

>^ ■ pulsars in the Galaxy above this luminosity limit. Using a pulsar current analysis, we 

J-j | derive the birth rate of this population to be 1.4±0.2 pulsars per century. An important 

conclusion from our work is that the inferred radial density function of pulsars depends 
strongly on the assumed distribution of free electrons in the Galaxy. As a result, any 
analyses using the most recent electron model of Cordes & Lazio predict a dearth 
of pulsars in the inner Galaxy. We show that this model can also bias the inferred 
pulsar scale height with respect to the Galactic plane. Combining our results with 
other Parkes multibeam surveys we find that the population is best described by an 
exponential distribution with a scale height of 330 pc. Surveys underway at Parkes 
and Arecibo are expected to improve the knowledge of the radial distribution outside 
the solar circle, and to discover several hundred new pulsars in the inner Galaxy. 

Key words: pulsars: general — pulsars: searches — pulsars: timing — stars: neutron 
— methods: statistical 



1 INTRODUCTION 

The Parkes multibeam pulsar survey of the Galactic plane is 
the most successful large scale search for pulsars so far un- 
dertaken. Five previous papers in this series have presented 
timing parameters for 600 newly discovered pulsars and have 
discussed various aspects of the survey results. A detailed 
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description of the 13-beam 20-cm receiver system and dat a 
reduction software can be found in Manches ter et alJ i200ll) . 
Prelimi nary remarks about the population stati stics were 
made bv lMorris et alJ l)2002l) . lKramer et alJ l|2003ft discussed 
the association of young pulsars with unidentified EGRET 
sources and gave evidenc e for pulsars tracing the spiral-arm 
structure of the Galaxy. iHobbs et, al.l l)2004Fl discussed the 
detection of previousl y known pulsars in the survey. Finally, 
iFaulkner et alJ (I2004T) implemented and discussed improved 
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processing and candidate selection strategies for the entire 
survey database with an emphasis on finding binary and 
millisecond pulsars. 

In this paper, we present the discovery parameters and 
provide timing solutions for a further 142 new pulsars. This 
brings the total number of pulsars found in the survey to 
742. Although future searches of the data may produce a 
few additional discoveries, the vast majority of "normal" 
(i.e. non-recycled) pulsars detectable by the survey have now 
been found. It is therefore appropriate to use this sample to 
place new constraints on the Galactic distribution and birth 
rate of normal pulsars. Of particular interest is the distribu- 
tion in Galactocentric radius. Due to severe selection effects 
on low-freque ncy (< 1 GHz) pulsar surv eys, earlier stud- 
ies (see, e.g.. iTavlor fc Manchester! Il977fl were somewhat 
hampered by small number statistics in the inner Galaxy 
and were often limite d to local population analyses (see, 
e.g., iLvne et al.lll99Sl) . In recent years, with the discovery 
of a larger number of pulsars in 1.4- GHz surveys there has 
been a growing body of evidence which suggests a deficit 
of pulsars in the inner Galaxy relative to a simple model 
wher e the radial density profile follows a Gaussian distribu - 
tion Jjohnstonll990l : lYusifov fc Kucukl2004t lLorimeill2004l) . 
One of the main goals of this paper is to use the new sample, 
which provides greatly improved pulsar statistics in the in- 
ner Galaxy, to provide an updated analysis on the Galactic 
distribution. 

The plan for this paper is as follows. In J^Jwe present 
the basic timing parameters, pulse widths, mean profiles and 
flux densities for the 142 new pulsars. The sample contains 
a number of interesting individual objects, including several 
binary and millisecond pulsars which are discussed in In 
521 we describe our Galactic population analysis. The main 
conclusions from this work are summarized in J5] 



2 DISCOVERY AND TIMING OF 142 
PULSARS 

The pulsars presented here were discovere d primarily us- 
ing th e processing schemes described by iFaulkner et alJ 
(2004). In brief, to compensate for the dispersive effects 
of the interstellar medium, data from each telescope beam 
were de-dispersed at 325 different trial values of disper- 
sion measure (DM) spanning the range to 2203 cm" 3 pc. 
The de-dispersed time series were then subjected in turn 
to: a standard fast Fourie r transform (FFT) ana lysis to 
search for periodic signals iManchester et al"ll200 lT) . a seg- 
men ted FFT search for ac celerated signals from binary pul- 
sars iFaulkner et all2004l) . a fast-folding algorithm (Kramer 
et al., in preparation) primarily sensitive to periods in 
the range 1-10 s and a search for dispersed single pulses 
jMcLauehlin et alJ [2006) . Pulsar candidates from all but 
the latter search were selected using a graphical tool which 
displayed various aspects of the search parameter space, 
e.g. signal-to-noise ratio (S/N) versus pulse period. Promis- 
ing candidates were noted for later re-observation at Parkes. 

Following the confir mation and position al refinement 
procedures described by iMorris et al] ll2002h . each pulsar 
was observed regularly for at least one year using one or 
more of the Parkes, Lovell and Arecibo telescopes. The 
observing systems used at Parkes and at Jodrell Bank 



are de scribed bv iManchester et al] (1200 1]) and lMorris et all 
(2002) respectively . The Arecibo timing is described by 
iHobbs et ail (120041) and IStairs et al] <200Sl) . For each pul- 
sar, pulse times of arrival (TOAs) were determined from 
the individual obse rvations using standard pulsar timing 
techniques (see, e.g.. lLorimer fc Kramerll2005f) impleme nted 
in the psrchive software package llH^tonet al]l2004|i 1 . A 
model containing the spin, astrometric and (if necessary) 
any binary parameters was fitted to the TOAs using the 
tempo timing package 2 . 

The positional information from these fits is provided in 
equatorial and Galactic coordinates in Tabled Subsequent 
columns in this table contain information on the discovery 
of each pulsar: the beam number (corresponding to the 13 
beams of the multibeam receiver) for the strongest discov- 
ery observation, the radial distance between the centre of 
this beam and the position of the pulsar (distances greater 
than one beam radius can occur if the pulsar scintillates or 
nulls or the closest pointing was contaminated by interfer- 
ence) and the S/N of the profile 3 during this observation. 
The observations used to form TOAs were added together 
to provide a characteristic pulse profile for each pulsar at 
1400 MHz (Figure [TJ. The final three columns in Table □ 
contain the flux densities measured from these mea n pro- 
files (for further details, see IManchester et al.ll200ll) . pulse 
widths at 50% and 10% of the pulse height. The 10% width 
is not measurable for pulsars with mean profiles that have 
poor S/N. For profiles containing multiple components the 
widths are measured across the entire profile. 

The rotational parameters from the timing analyses 
are given in Table |2] In column order, this table provides 
each pulsar's name, Solar-system barycentric pulse period, 
period derivative, epoch of pulse period, the number of 
TOAs used in the timing solution, the MJD range covered 
by the timing observations, root-mean-square value of the 
observed— model residuals and the DM. The data have been 
folded in turn at two and three times the tabulated periods 
to confirm that they represent the fundamental periods of 
the pulsars, rather than a harmonic. 

Various derived parameters for the new pulsars are pre- 
sented in Table |3] For each pulsar, we list the base-10 log- 
arithms of the characteristic age, surface dipole magnetic 
field strength and rate of loss of rotational energy. The fi- 
nal columns contain the pulsar distances, height above the 
Galactic plane and luminosities 4 . The distances are com- 
putcd from their DMs assuming the lTavlor fc Cordesl (Il993l) 
and ICordes fc Laziol (1200211 models for the Galactic distri- 
bution of free electrons. The former is used for consistency 
with earlier papers in this series. The latter is more up-to- 
date and used for the population analysis in H4.2I 



1 http : //psrchive . sourcef orge .net 

2 http : //www. atnf . csiro . au/research/pulsar /tempo 

3 Throughout this paper we discuss S/N measurements made 
in the time domain, i.e. determined directly from the in tegrated 
pulse profile (for further details, see lLorimer fc Kramej|2005l) 

4 Note that throughout this paper we will use the term luminosity 
L to describe the quantity L = Sd 2 , where S is the flux density 
and d is the distance. As this omits any geometrical and beaming 
factors, t his quantity is often refer red to as a pseudoluminosity 
(see, e.g,. lArzoumanian et al1l2002l) . 
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3 DISCUSSION OF INDIVIDUAL OBJECTS 

The 7.7-s pulsar J1001— 5939 has the longest spin period 
in our sample, and is the second longest kno wn for a ra- 
dio p ulsar after the 8.5-s pulsar J2144— 3933 ijYoung et alJ 
1999). Unlike PSR J2144-3933, the much larger period 
derivative measured for J1001— 5939 means that the pul- 
sar lies above the "death lin e" in the P — P diagram 
iRuderman fc Sutherlandll975l) . However, its implied dipole 
surface magnetic field is below the quantum critical limit 
jBaring fc Hardind liool . As a result, PSR J1001-5939 
does not pose any serious problems to theories of radio emis- 
sion. Although the initial survey detection of this pulsar 
was registered in its second harmonic in the standard FFT 
search, its detection was much more significant using the 
fast-folding algorithm. Full details of this aspect of the sur- 
vey data processing will appear elsewhere (Kramer et al., in 
preparation). 

In addition to the new discoveries of rotating ra- 
dio transient sources i n the single-pulse search described 
by iMcLaughlin et alJ |2006), one further pulsar, PSR 
J1624— 4613, was only detectable in this part of the anal- 
ysis. Our follow-up observations show that this pulsar is 
only detectable about 70% of the time. This and the lack 
of detection by the periodicity analysis suggests that this 
pulsar spends a significant fraction of time in a null state. 
The discovery of such pulsars alongside the radio transients 
highlights the value of the single-pulse search analysis as an 
important component of pulsar survey pipelines. 

The youngest pulsar in our sampl e, PSR J1357— 6429, 
has been extensively discussed bv lCamilo et al.N2004h . This 
pulsar, which has suffered a Vela-like glitch, is possibly as- 
sociated with the radio supernova remnant G309.8— 2.6. We 
have searched for coincidences between the newly discov- 
ered pulsars that lie within the bounds of the 231 supernova 
remnants in Green's catalogue as well as the 35 rem nants 
recently found in the inner Galaxy jBrogan et alj feooeV Al- 
though five additional matches were found (G4.2— 3.5 and 
J1808-2701; G27.8+0.6 and J1839-0436; G321.9-0.3 and 
J1519-5734; G343. 0-6.0 and J1718-4539; and G343. 0-6.0 
J 1724— 4500) none of the pulsars have characteristic ages 
below 6 x 10 5 yr or an angular offset less than 0.45 rem- 
nant radii. Such apparent coincidences are common and ex- 
pected, given the high density of pulsars and supernova rem- 
nants along the Gal actic plane f Ga ensler fc Johnston 1995; 
lLorimer et alJ [l998) . To demonstrate this, we artificially 
shifted the remnant centres by 1° in Galactic longitude. 
The subsequent fake cross-correlation analysis resulted in 6 
matches with similarly old pulsars and large angular offsets. 
We therefore do not believe any of the above coincidences 
are genuine pulsar-supernova remnant associations. 

Several of the profiles shown in Fig. exhibit the char- 
acteristic hallmark of interstellar scattering, i.e. a one-sided 
exponential tail. The most extreme example in the cur- 
rent sample is PSR J1901+0435, a 691-ms pulsar located 
at Galactic longitude / = 38°, which shows a scattering tail 
encompassing most of its rotational phase. Furt her study of 
this and other pulsars in the sample (see, e.g., iBhat et alJ 
120041) will help improve models of the distribution of scat- 
tering material in the interstellar medium. These detections 
of highly scattered pulsars also suggest that, even at the 
survey frequency of 1.4 GHz, a number of nominally bright 



pulsars are not detectable due to scattering. Further surveys 
of the Galactic plane at even higher observing frequencies 
are required to detect these pulsars. 

Our sample includes the two solitary millisecond pulsars 
J1 801-1417 and J 1911+ 1347, both originally announced 
by Faulkn er et alJ (j?004). Improved timing parameters for 
these pulsars are given in Tables and [5] These discoveries 
bring the sample of solitary millisecond pulsars known in the 
Galactic disk to 15. In keeping with the other isolated mil- 
lisecond pulsars, these new sources have low radio luminosi- 
ties compared to their counterparts with binary companions 
jBailes_et_aljJr997t jKramer et al1ll998l) . As noted recently 
bv lLommeriet^ilT(l2006l) the difference in scale heights of the 
binary and isolated pulsars, despite their having statistically 
identical velocity dispersions, also points to a difference in 
luminosities for the two populations. It is currently unclear 
why the luminosity functions differ for the two populations. 

Our sample includes 13 binary pulsars, the properties 
of which are summarized in Table 01 Of these, seven are 
discussed elsewhere: PSR J1638— 4725 is a long-period bi- 
nary pulsar in a highly eccentric orbit around a massive 
companion (Lyne et al., in preparation), PSR J1744— 3922 
is a short-period bina ry system with a low-mass companion 
jFaulkner et al.l2004l). PSR J1802-2124 has a companion of 
intermediate mass jFaulkner et alll2004h. PSR J1756- 2251 
is a double-neutron-star binary jFaulkner et al.ll2005Tl and 
PSRs J1751-2857, J1853+1303 and J 1910+1256 are mi l- 
lisecond pulsars in wide binary orbits JStairs et al.ll2005ri . 
In Tables Qj] and |S| we provide orbital parameters for the re- 
maining six binary systems. 

PSRs J1711-4322 and J1822-0848 listed in Table |U 
appear to be similar in character to PSR B0820+02; i.e. a 
long-period pulsar in a long-period orbit about a low-mass 
white dwarf companion. PSRs J1125-6014, J1216-6410 
and J1841+0130 are typical of other recycled pulsar bi- 
nary systems with orbital periods of typically a few days 
and low-mass white dwarf companions. All of these bi- 
naries follow the spin period versus orbital period and 
ma gnetic field versus orbital period c orrelations discussed 
by Ivan den Heuvel fc Bitzarakil 1119951) which support the 
largely empirical hypothesis that accretion onto the neutron 
star causes a reduction in the magnetic fie ld strength (see, 
e.g.. IBhattacharva fc van den Heuvellll99ll) . 

Another important tracer for the evolution of binary 
pulsars is the orbital period versus eccentricity co rrelation. 
Accor ding to the fluctuation-dissipation theory of IPhinnevI 
(1992), a relationship between these parameters is expected 
for all pulsar-white dwarf binary systems which have under- 
gone an extensive period of mass transfer. With the excep- 
tion of PSRs J1841+0130 and J1853+1303, which deviate 
from the theoretical prediction by over an order of magni- 
tude, the pulsar-white dwarf syst ems in our sample follow 
the P;, — e correlation. As noted bv lStairs et al.l (120051) these 
two exceptions may simply be indicative of a natural scatter 
about the relationship. Further work is required to establish 
whether these deviations are indicative of a different evolu- 
tionary channel, such as the formation of the neutron star 
via accretion induced collap se of a white dwarf in an ac cret- 
ing binary system (see, e.g.. INaravan fc Pophamlll989t) . 

Slightly less typical is the 28-ms pulsar, PSR 
J1439— 5501. This binary pulsar with a circular orbit 
of period 2 days has a minimum companion mass of 
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1.1 Mq. It appears to be among the growing class of re- 
cycled pulsars with more massive CO white dwarf com- 
panions, the so-call ed "intermediate-mass binary pulsars" 
icamilo et al"lll996ft . Its closest known counterpart among 
the known systems is th e 43-ms binary pulsar J1157— 5112 
(Edw ards fc Ba ilcs 2001). Another likely intermediate-mass 
system is J1802 — 2124 for which the minimum companion 
mass is 0.8 M (iFaulkner et al.ll2004l) . 

Related to the binary pulsars is the isolated pulsar 
PSR J1753— 1914. With a spin period of 63 ms and a spin- 
down rate of only 2 x 10 -18 , this pulsar appears to be- 
long to the class of so-called "disrupted recycled pulsars" 
(lLorimer et alJl2004h . Since these objects have spin proper- 
ties similar to the recycled pulsars in double-neu t ron-st ar bi- 
nary systems, it was proposed bv lCamilo et al.l (Il993l) that 
they are the result of the systems which di srupted at the 
time o f the second supernova. As discussed bv lLorimer et all 
(2004), there appears to be a deficit of these systems among 
the known pulsar sample. Further investigations, especially 
using binary population syntheses, are required to investi- 
gate the frequency of these isolated recycled pulsars relative 
to double neutron star systems. 



4 THE GALACTIC PULSAR POPULATION 

To date, the total number of pulsars discovered by the 
Parkes multibeam (hereafter PMB) survey is 742. This sub- 
stantial haul, together with a further 11 transient radio 
sources found in the survey data asso ciated with rotating 
neutron stars (iMcLauehlin et al . 2006), means that our sur- 
vey is a phenomenal new probe of the Galactic neutron-star 
population. We also make use of the recently published high- 
latitude pulsar survey of the regi on —140° < I < —10 0° and 
|6| < 60° (hereafter PH survey. iBurgav et aljl2006T) . Since 
the PMB and PH surveys were carried out using the same 
observing system, the pulsars detected represent a reason- 
ably homogeneous sample from which to proceed. 

Although the total number of pulsars detected by the 
PMB and PH surveys is 1055, for this analysis, where we 
are concerned with the properties of normal pulsars, we use 
only those 1008 non-binary pulsars with P > 30 ms. In this 
sample, 976 pulsars were detected solely by the PMB sur- 
vey, 22 were detected solely by the PH survey and 12 were 
detected jointly. While there are two additi onal large-scale 
surveys using the P arkes multibeam system llEdwards et all 
bOOltlJacobvl BoOS). these are not included at this stage of 
the analysis due to lack of flux densities for many of the pul- 
sars from these surveys. However, these surveys are utilized 
later in the discussion ( H4.5.2H when the scale height of the 
population is investigated. 

In the following, we investigate various aspects of the 
Galactic pulsar population after accounting, as far as possi- 
ble, for the known selection effects. We begin by presenting 
a model for the survey detection process which is in good 
agreement with the observed S /Ns from the PMB survey. We 
then use this detection model in a Monte Carlo simulation 
to iteratively deduce various underlying Galactic population 
properties. Finally, we apply a pulsar current analysis on our 
optimal model to deduce the birth rate of the population. 



4.1 Survey detection model 

Although iManchester et alj (l200lT) presented a detailed for- 
malism to calculate the survey sensitivity, we have found 
that a si mpler approach is t o use the radiometer equation 
(see, e.g., iDewev et alJ ll984) . As we show below this pro- 
vides a self-consistent model for t he detection statistics. This 
appro ach was also adopted by iFaucher-Giguere fc KaspH 
( 2006) where it was noted that there app ears to be some in- 
consis tencies with the procedure given bv IManchester et alJ 
(2001). The radiometer equation gives the observed S/N in 
terms of various pulsar and system parameters as follows: 

Jn v Avr I P -W 
S/N^ 14 oo G ^^ — . (1) 

Here S1400 is the pulsar mean flux density at 1400 MHz 
(mjy), G is the effective telescope gain (K/Jy), n v is the 
number of polarizations summed, Av is the observing band- 
width (MHz), r is the integration time (s), (3 accounts for 
S/N losses, T is the system temperature (K), P is the pulse 
period (s) and W is the observed pulse width (s). Since both 
the PMB and PH surveys used the same telescope, receiver 
and data acquisition system, the only parameter that is dif- 
ferent is the integration time r = 2097 s for the PMB survey 
and 262 s for the PH survey. The con stant parameters are 
n v = 2, Au = 288 MHz and j3 = y^T^ c ± 1.25 for S/N losses 
due to one bit digitization (see, e.g., lLorimer fc Kramer! 
2005). Gain and system temperature values depend primar- 
ily on which element of the 13-beam receiver system was 
used to make the dete ction. We use the values given by 
IManchester et all (1200 if) . Sky temperature values vary as a 
function of Gal actic latitude and lon gitude. To account for 
this, we use the lHaslam et ail (Il982f) 430-MHz all-sky cata- 
logue scaled to the observing frequency of 140 MHz assum- 
ing a^ sky background spectral index of —2.6 dLawson et alJ 

To account for the drop in flux density away from the 
beam centre, we ca n assume a Gaussian beam (see, e.g., 
lLorimer et"al]ll993ft to write 

G = Go exp(-2.77r 2 /^ 2 ), (2) 

where Go is the telescope gain at the beam centre, r is 
the offset from the beam centre and w is the full-width at 
half power of the telescope beam. For each detection un- 
der consideration, we ad opt the Go and w values given by 
IManchester et aTl ||2001) for the appropriate beam of the 
multibeam receiver when making this correction. To calcu- 
late the equivalent pulse width W from the observed 50% 
width, W50, we make a simplification that the pulse shape is 
Gaussian in form (a reasonable approximation for most pul- 
sars). For a Gaussian pulse with an intensity of unity and a 
standard deviation a, we find that W50 = cta/8 In 2 and W is 
just the area under the pulse (i.e. W = ayphi'). Eliminating 
a from these two equalities yields the relationship 

W = Wso^/j^2 a 1.06W 50 . (3) 

Figure |21 shows a comparison between the modeled and ob- 
served S/N values using data taken from Tableland the 
corresponding data from the earlier papers in this series. 
The error bars are derived directly from the uncertainties in 
the measured flux densities of each pulsar. As can be seen, 
the agreement between our model and observed S/N values 
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is good, being largely free from systematic trends and with 
a scatter that is not much larger than the uncertainties. In 
addition to providing a sanity check on the system perfor- 
mance, this analysis demonstrates that Equation provides 
a good model of the survey detection thresholds and, for the 
sample of normal pulsars presented here, we recommend its 
use in future analyses. 



4.2 Galactic distribution determination method 

Using the above model for the PMB and PH surveys, we 
now describe a method to derive the underlying probability 
density functions (PDFs) in pulse period (P), 1400-MHz 
radio luminosity (L), Galactocentric radius (R) and height 
above the Galactic plane (z). The modeling techni que we 
use is based on the iterative proced ure developed by jLargd 
Jl97ll) . lTavlor fc Manchester! dl977l) and lLvne et"alTil985t) . 
A preliminary acc ount of the analysis below was presented 
bv lLorimerl J2004D . 

Our goal is to produce a model of the Galactic pulsar 
population which, when filtered through the above survey 
detection models of the PMB and PH surveys produces a 
set of "model detectable pulsars" which closely matches the 
observed sample. A Monte Carlo simulation package 5 was 
developed for this work and is freely available for checking 
the results presented here, and for further use. 

We simulate the Galactic pulsar distribution by draw- 
ing pulsars from independent numerical PDFs in R, z, L and 
P. The justification for assuming no correlation between the 
distributions is twofold. Firstly, no statistically significant 
correlati ons between t hese parameters exist in the known 
sample iLorimeJ l2004Tl . Secondly, our modeling philosophy 
is to create a snapshot of the currently observable pulsar- 
population, rath er than attempting to follow th e time evo- 
lution (see, e.g., IFaucher-Gieuere fc KaspHl2006t and refer- 
ences therein) . Although there are advantages to taking the 
time- evolution approach, particularly when investigating the 
evolution of pulsar luminosities and magnetic field strengths, 
we prefer to adopt the simpler procedure here which concen- 
trates on issues relating to the present-day spatial, period 
and luminosity distributions. 

We choose the four numerical PDFs as follows: for the 
R distribution we consider 15 equal zones between and 
15 kpc; for the z distribution we consider 31 equal zones be- 
tween -1.5 and 1.5 kpc; for the L distribution we consider 12 
logarithmically spaced zones between 0.1 and 1000 mjy kpc 2 
each one third of a decade wide; for the P distribution we 
consider 15 logarithmically spaced zones between 30 ms and 
10 s each one sixth of a decade wide. To begin with, we as- 
sume no prior knowledge and distribute pulsars completely 
uniformly in space and use flat PDFs in log L and log P. We 
also experimented with other starting PDFs and concluded 
that any function which reasonably samples the parameter 
range produces results that are consistent with those found 
below. 

During the course of this work, we found that the most 
important factor limiting our determination of the true spa- 
tial density of pulsars in the Galaxy was the uncertainty 
in the Galactic distribution of free electrons. We therefore 
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distinguish between two different models. The first model 
(hereafter referred to as "S" for smooth), treats the pulsar 
and free electron density distributions as smooth and az- 
imuthally symmetric functions with respect to the Galactic 
centre. The second model (hereafter referred to as "C" for 
clumpy), attempts to incorporate Galactic spiral arm struc- 
ture and a heterogeneous distribution of free electrons. 

To specify the location of any model pulsar, we use 
a Cartesian (x, y, z) coordinate system where the Galactic 
centre is at the origin, and the Sun is at (0.0,8.5,0.0) kpc. 
With this definition R = \/ x 2 + y 2 . For a given pair of R 
and z values, then, we can calculate the Cartesian coordi- 
nates of each model pulsar. We achieve this in two differ- 
ent ways. In model S, we choose a random azimuthal angle 
in the disk of the galaxy 8 from a flat PDF in the range 
< < 2ir. The coordinates are then simply x — RsinO 
and y = R cosO. In model C, we follow the procedure de- 
scribed by IFaucher-Gieuere fc KaspH (j200G) to weight the 
pulsar distribution along model spiral arms and calculate 
the appropriate x and y positions. In both cases, after de- 
termining the x, y and z coordinates of each pulsar, it is 
straightforward to compute its true distance from the Sun 
d and apparent Galactic coordinates, I and b. Finally from 
the distance and luminosity of each pulsar, and the definition 
L — Snood 2 , we calculate the true 1400-MHz flux density, 

Sl400- 

As mentioned in the previous section, the true flux den- 
sity needs to be corrected for the position offset between 
a pulsar and the nearest telescope pointing in the survey. 
The PMB survey comprises data from 40077 beam positions 
shown in Figure [3] Although the pointings principally cover 
the longitude range —100° < I < 50° and |6| < 5°, as can be 
seen there is some coverage of the region —140° < I < —100° 
and |fc| < 3°. From the sky position of each model pulsar we 
use this database to find the closest beam position in the sur- 
vey. To model the PH survey, which had a total of 6456 beam 
positions in the region -140° < I < -100° and |6| < 60°, we 
randomly choose an offset within the full-width half power 
point of a telescope beam. From the corresponding position 
offset and true flux density, we use Equation to find the 
effective flux density, S. 

To calculate DMs, we integrate two different models 
for the free electron density out along the appropriate line 
of sight defined by I and b until the distance d is reached. 
For model S we u se the azimuthally symmetric and smooth 
model 6 derived bv lLvne et alJ lll985l) . while for model C, we 
adopt the free electron distribution due to lCordes fc Laziol 
( 2002) which attempts to account for the clumpy nature of 
the interstel lar medium includ ing s piral arm structure. A s 
discussed bv lLvne et alJ (119851) and lCordes fc Laziol J200^ . 
the overall random errors in both models result in distance 
uncertainties of the order of 20%. We account for this in 
both models by dithering all true distances using a Gaus- 
sian distribution with a fractional error of 20% and use the 
symbol * to denote parameters affected by this distance un- 
certainty. Hence, in additi on to th e dithered distance, d, 
we have L = Sd 2 , R = \J x 2 + y 2 and z — ds'mb as the 



° We note that a m ore recent electron density model proposed by 
iGomez et alJ was not used here due to its lack of predictive 

power for pulsars in the inner Galaxy (R < 4 kpc). 
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corresponding radius, R, and z-height estimates. For self- 
consistency, when comparing our model samples with the 
observations below, we use the model L, R and z quantities. 

Having computed the position and flux density of each 
pulsar, we are now in a position to model its detectability. 
To determine the model S/N using Equation Q we need to 
find the detected pulse width, W. This can be written in 
terms of the intrinsic pulse width Wxat and a number of 
other factors in the following quadrature sum: 



w = VwL + wEm + ^s. 



(4) 



Here Wdm is the dispersive broadening across individual fil- 
terbank channels, t sa m P is the data sampling interval (mul- 
tiples of 250 fis for the PMB survey and multiples of 125 
jj,s for the PH survey) and iscatt is the interstellar scattering 
time. In the following, we examine each of the quantities on 
the right hand side of this equat ion in turn. 

It is well known (see, e.g., iLvne fc Manchester! ll98sF ) 
that pulse widths are correlated with pulse period. During 
the course of this work, it became obvious that the observed 
scatter in the W-P plane (see Fig. cannot be explained 
by a random relationship. Some correlation between W and 
P is to be expected, given the known var iation in beam 
size with pulse period (see, e.g.. lBiggslll990f) . despite being 
blurred by an arbitrary viewing geometry for each pulsar. 
Rather than adopt a complicated beam model, we adopt a 
purely empirical approach in our simulations by assigning 
the intrinsic pulse width, Wint, as follows: 



log Wint = log 



0.06 



(-y 

Vms/ 



(5) 



where T is a suitably normalized Gaussian distribution 
about the origin with a standard deviation of 0.3. These 
parameters were found by experimentation so that, after 
taking into account the various additional pulse broadening 
effects described below, detected pulse widths of the model 
pulsars are in reasonable agreement with the observations 
(see Fig. 2J. Although the exact values of the above nu- 
merical factors are not crucial to our final results, models 
without this correlation built in were found to produce very 
unsatisfactory comparisons with the real sample. 

Turning now to the instrumental terms in Equation 01 
we can insert the parameters of both surveys and write the 
dispersive broadening across individual 3-MHz frequency 
channels as follows: 



Wdm = 9.6 /xs 



DM 



' pc 



(6) 



As described by M anchester et ai] i200lf) . the sam- 
pling time tsamp depends on DM tr iai- A file containing 
a list of trial DMs and sampling times is available at 
http:/ /www. blackwellpublishing.com/products/journals/su 
ppmat /MNR/.... /PMSURV.DMs 

Finally, in order to model interstellar scattering, we 
make use of the well-known correla tion between DM and the 
scattering timescale t SC att (see, e.g.. lSlee et alJl98o[). Scaling 
the be st-fit relationship from the recent study of iBhat et all 
J2004D to 1400 MHz, we find 

(7) 



T = 0.154D + 1.07V - 7, 



where T and T> are respectively the base- 10 logarithms of 
tscatt (in ms) and DM (in cm -3 pc). To account for the 



dispersion about this relationship (see Fig. 4 of Bhat et 
al. 2004), we dither the model scattering times using a nor- 
mal distribution around T with a standard deviation of 0.8. 

The above discussion has detailed all the main steps 
in the simulation which allow us to calculate a predicted 
S/N for each model pulsar using Equation Those pulsars 
with S/Ns greater than 9 were deemed to be theoretically 
detectable. Each simulation proceeded until 25,000 model 
detectable pulsars were accumulated. This number was cho- 
sen to be much larger than the t rue sample size to minimiz e 
statistical fluctuations (see, e.g.. iBhattacharva et aljll992Tl . 
For each simulation, histograms of the observed distribu- 
tions of P, L, R and z were formed and normalized to have 
the same area as the sample of real observed pulsars. 

To assess the overall quality of a simulation, we calcu- 
lated the reduced \ 2 statistic summed over the bins of the 
P, L, R and z distributions. We note however that \ 2 is 
used here only as a figure-of-merit, and not as a means for 
hypothesis testing or rejection. For a typical starting model 
with flat PDFs in R, L, z and P the reduced \ 2 was, as 
expected, rather large, with typical values of order 1000. To 
improve these models, for each distribution, we computed a 
set of correction factors with one factor per bin defined as 



Ci 



Ri-Mj 
Mi 



(8) 



where Ri and M; are respectively the real and model ob- 
served number of pulsars in the i th bin of one of the P, L, 
R and z histograms. A new set of input distributions was 
then generated by multiplying each bin by the correction 
factor and adding this product to the original bin value. For 
example, the new z distribution, z , is computed as follows: 



Zi -\- Zi X Ci 



(9) 



The resulting new set of PDFs should provide a better match 
to the observations. The simulation is then re-run and a new 
set of correction factors computed. We found that this pro- 
cedure rapidly converged to produce a model with a much 
lower reduced x 2 (typically of order unity) that stabilized af- 
ter only a few iterations. We note in passing here that there 
appears to be very little covariance between the four PDFs. 
We found, through experimentation, that it was not possi- 
ble to reach any satisfactory convergence by keeping one or 
more parameters constant during the iterations and varying 
the others. In this regard, at least within the overall limita- 
tions of the method we are using, we are confident that the 
models presented below are optimal and do not represent 
'local minima' solutions. 



4.3 A self-consistency test of the model 

To verify that the above procedure can successfully re- 
cover the various Galactic distribution properties, we sim- 
ulated two fake populations with different radial density 
functions. In the first population, pul sars were drawn from 
a simple Gaussian profile (see, e.g.. iNaravanl IT987fl with 
a standard deviation of 6.5 kpc. For the second pop- 
ulation we adop t ed th e more recent profile derived by 
lYusifov fc Kuciikl (|2004) which peaks at ~ 4 kpc and has 
a deficit of pulsars in the inner Galaxy. In both simula- 
tions, we assumed an exponential z distribution with a scale 
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height of 350 pc jMdzinarishvili fc Melikidzd [2004 ) , a log- 4.5 Discussion of the results 



normal distribution of luminosities with a mean of -1.1 
and standard deviation 0.9 (b oth in log 10 of mjy kpc 2 , 
IFaucher-Gieuere fc K aspi 2006) and a gamma function for 
the spin periods which peaks at ~ 600 ms as described by 
iGil fc Hani lll996t) . The assumed pulse width versus period 
relationship was given in the previous section. Both models 
were run such that 1005 model pulsars were detected by our 
models of the PMB and PH surveys. The resulting model 
observed samples were then used as fake observed popula- 
tions, i.e. for the purposes of this test, they take the place 
of the real observed sample. 

After applying the Monte Carlo procedure described in 
the previous section to these fake samples, we were able to 
reproduce the functional form and correct number density 
over the parameter space of each model. Fig. shows the 
results in terms of the input and derived radial density func- 
tions. These tests give us confidence that, when applied to 
the real observed sample of pulsars, our procedure provides a 
reliable means of estimating the true spatial and luminosity 
functions of pulsars in the Galaxy. 



4.4 Application to the observed pulsar sample 

We now apply the above procedure to the real sample of 
1008 isolated pulsars detected by the PMB and PH sur- 
veys. As discussed above we consider two diff erent cases: 
a "sm ooth" (S) model of the Galaxy using the iLvne et alJ 
(1985) electron density distribution a nd a "clumpy" (C) 
model using the lCordes fc Laziol i2002f) electron density dis- 
tribution and following the spira l arm modeling described by 
IFaucher-Gieuere fc Kaspil i2006t) . In each case we computed 
distances to both the real and model pulsars using the ap- 
propriate electron density model. 

Figs.|S|and[7|show the real observed and resulting model 
underlying distribution functions derived from our analysis. 
The error bars shown here are purely statistical estimates 
with a fractional uncertainty equal to l/y/N, where N is 
the number of observed pulsars in each bin. To parameter- 
ize our results in a form that may be of use to others, we 
also attempt to fit various analytic functions to the under- 
lying distributions; these are shown by the smooth curves in 
each figure. For the radial density profile, we use a gamma 
function 



<7^ 



exp 



-C 



R — R, 



o 



R 







(10) 



for the z distribution, we use the exponential function 
N = Dexp(-\z\/E); (11) 
for the luminosity distribution, we fit a simple power law 
logN = FlogL + G; (12) 
finally, we fit the period distribution to the function 
[logP-/] 2 " 



AT (log P) = Hexp 



2J 2 



(13) 



Table[7|summarizes the resulting fit parameters A through J 
for both models. The dotted lines in the figures compare our 
results with various other distributions and are described in 
turn below. 



4-5.1 Radial distribution 

The underlying radial density distribution of each model 
provides a straightforward means of calculating the implied 
number of potentially detectable pulsars (i.e., above the lu- 
minosity threshold of 0.1 mjy kpc 2 and beaming toward 
us) in the Galaxy, Ng- For models S and C respectively we 
find N G = 28000 ± 1100 a nd 30000 ± 1100. Adopting the 
iTauris fc Manchester! (I1998T) beaming model, these numbers 
translate to 148000 ± 6000 and 155000 ± 6000 active ra- 
dio pulsars in the Galaxy with luminosities above 0.1 mjy 
kpc 2 . While these numbers are consistent, the shape of the 
radial density distributions are quite different. As seen in 
Figs. Hjji and^K, where the dotted lines show the radial dis- 
tribution of free electrons, we find that the pulsar distribu- 
tion very closely follows that of the electrons in both cases. 
This striking correlation can be understood when one real- 
izes that model C has an enhancement of electrons in the 
range 3 < R < 5 kpc. Thus, at low Galactic longitudes, the 
model will preferentially place pulsars either on the near 
or far side of the Galaxy, close to the electron density en- 
hancement. For model S, where the electron density peaks 
at R = 0, such pulsars are preferentially placed towards 
the centre of the Galaxy. These different radial distributions 
have implications for the detection rates of future surveys 
as we discuss later in H4.7I 

We conclude from the above discussion that the Galac- 
tic pulsar and electron density distributions in the inner 
Galaxy are highly covariant. Bas ed on our cur rent state of 
knowledge, i.e. that the lCordes fc Laziol (l2002f) model pro- 
vides a far better description of the free electron distribu- 
tion than all of its predecessors, our optimal model for the 
Galactic pulsar population is model C. However, given the 
uncertainty in the electron distribution in the inner Galaxy, 
care should be taken when interpreting the results of our fit 
which shows a peak in the pulsar distribution at 3-5 kpc. 
Further progress on the radial density distribution of pulsars 
requires independent distance estimates for more pulsars in 
the inner Galactic quadrants. Another possibility worth ex- 
ploring, though beyond the scope of this paper, is to treat 
the electron density distribution as a free parameter in the 
Monte Carlo simulations. For now, we simply caution that 
the deficit of pulsars in the inner Galaxy noted here (see 
also the earlier pap ers bv lJohnstonlll99'4l lYusifov fc Kuciikl 
120041 : lLorimerll2004l) depends on the assumption of a corre- 
sponding deficit in the free electron density. 

4-5.2 z distribution 

We now consider an interesting result from our fits to 
the z distributions of both models, where it is seen that 
model S results in a scale height of 330 pc compared to 
only 180 pc for model C. The latter value is much lower 
than the canonical value of 300-350 pc expected a-priori 
from an independent study of th e local pulsar population 
jMdzinarishvili fc MelikidzejF2 004 ) and shown by the dotted 
line in Figs. |Sp and |7p. As for the radial distribution, we 
believe this difference to be a direct result of the different 
electron density distributions used in the two models. From 
a comparison of Figs.|Sp and[7p, we see that the z heights of 
the observed pulsar sample computed using the lLvne et alJ 
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( 1985) elec tron density m o del ar e much broader than when 
using the lCordes fc Lazid J2002f) model. 

The discrepancy is further highlighted when one uti- 
lizes the two other surveys carried out with the Parkes 
multibeam system, i.e. the Swinburne intermediate and high 
latitude pulsar surveys (respective ly SIL and SHL surveys, 
lEdwards et"aHl200lt Ijacobvll2005l) . Along with the PH sur- 
vey, these Swinburne surveys sample the distribution away 
from the Galactic plane and are therefore very sensitive to 
the z scale height of the pulsar population. As mentioned 
previously, these surveys were not included in the main anal- 
ysis because the majority of pulsars from them do not have 
measured flux densities. However, our models should be able 
to reproduce the relative yield of the PMB and PH surveys 
alongside the SIL and SHL surveys. We have carried out 
such an analysis and compare the observed pulsar samples 
from these surveys with the various models in Tabled] As 
can be seen, the expected survey detections for model S are 
generally in good agreement with the actual survey results. 
For model C, however, the predicted yield for all surveys 
away from the Galactic plane is much lower than observed. 

As a final sanity check, we have verified through ad- 
ditional simulations similar to those described in iJ4.3l that 
our modeling procedure is capable of reproducing any rea- 
sonable input scale height distribution. Based on these re- 
sults, we concl ude that there i s a po tential inconsistency 
when using the lCordes fc Lazid (|2002) model to derive the 
scale height of pulsars and consider a more likely value to 
be that found by model S, i.e. 330 pc. The apparent clus- 
tering of the PMB pulsars towa rds small z heights when 
using the ICordes fc Lazid l l2002h model was also noted by 
iKramer et al J J2003tl . It appears that further investigation of 
the z-distribut ion in the electron density model i s now war- 
r ant ed. While iMdzinarishvili fc Melikidzel ((2004) also used 
the ICordes fc~LarioT " (l2"o"o"2Tr electron density model, their 
analysis only considered pulsars within a cylinder of radius 
3 kpc centred on the Sun. This local population analysis 
was done to minimize selection effects. Although our analy- 
sis attemp ts to correct for these ef fects as far as possible, we 
believe the lCordes fc Lazid J2002h distance model is biasing 
our results towards an artificially lower scale height. 

In an attempt to correct for this bias, we have intro- 
duced a modified version of model C, denoted by model C', 
which is identical with the original model with the exception 
of having an exponential z distribution with a scale height of 
330 pc. As can be seen from the yields of this model in Table 
|3] the revised scale height is in much better agreement with 
pulsar surveys carried out at higher Galactic latitudes. We 
return to this model for the discussions on pulsar current 
analysis and survey predictions (sections 14.61 and 14.71 . 



4-5.3 Luminosity distribution 

For both models S and C, we find that the simple power 
law given in Equation 1121 provides an adequate description 
of the underlying luminosity function for L > 0.1 mjy kpc 2 . 
However, the slope of the distribution (-0.6 for model S 
and -0.8 for model C) is somewhat flatter than the canon- 
ical value d log Nj d log L = — 1 derived in early studies o f 
the pulsar population (see, e.g.. lTavlor fc ManchesteJl977l) . 
Overlaid on both Figs. [(J and [7J: is the log-normal fit 
to the luminosity distribution from the optimal model of 



iFaucher-Gieuere fc Kaspil (I2006I) . As can be seen, given the 
uncertainties involved, both mod els S and C are in reason- 
able a greement with the results of lFaucher-Giguere fc Kaspil 
(2006). 

The behaviour of the luminosity function below 
0.1 mjy kpc 2 is not well constrained by our analysis. Al- 
though we began our study using a luminosity distribution 
which extended to 0.01 mjy kpc 2 , we found that the few pul- 
sars in the observed sample with L < 0.1 mjy kpc 2 could be 
accounted for simply by a luminosity function with a lower 
limit of 0.1 mjy kpc 2 and the dithering of luminosities re- 
sulting from a 20% distance uncertainty. Better constraints 
on the shape of the low end of the pulsar luminosity func- 
tion should be possible from an analysis of deep targeted 
searc hes for young pulsars (for a recent review, see ICamilcl 
2004). 



4-5.4 Period distribution 

In Figs.|SJi and|7Jl, we compare our p eriod distribution with 
the parent distribution favoured by iKolonko et alJ (120041) 
in their analysis of pulse-width statistics. As can be seen, 
our distribution is significantly different, with a much higher 
proportion of short-period pulsars. Given the effort to quan- 
tify and account for period-dependent selection effects in this 
paper, we believe our fit using the log-normal distribution 
given in Equation 1131 provides a much better description of 
the parent population of normal pulsars and recommend its 
adoption in future pulse- width analyses. 

4.6 Pulsar current analysis 

To compute the birth rat e of the population, w e fol- 

Il98ll: 



Vivekanand & Naravanl |l9^|: 


Lorimer et alJ 


Lvne et al.l Il998l: IVranesevic et al 


2004) and 



1993 



a pulsar current ana lysis. We closely follow the approach 
described in detail bv lvranesevic et alJ (120041) and calculate 
a weight or "scale factor" for each of the 1008 pulsars in our 
main sample. The analysis proceeds by placing each of these 
pulsars in every position of our model galaxy and recording 
the number of detections made by our model PMB and PH 
surveys. The model galaxy used for these calculations was 
model C', i.e. the modified version of model C which has 
a z scale height of 330 pc and produces the best match to 
the other survey detection rates (see ^4.5.21 . To compute 
the apparent flux density, we use the observed luminosity 
of each pulsar and apply the empirical relationship given 
in Equation to calculate an intrinsic pulse width at each 
position in the model Galaxy. The scale factor is then the 
ratio of the total number of positions tried to the number 
of detections found. 

Having obtained a set of scale factors, the current or 
flow of pulsars from short to longer periods can be calcu- 
lated by binning the sample into period intervals of width 
AP. For a given period bin con taining n psr pulsars, following 

we define the pulsar current 



(14) 



IVivekanand fc Naravanl jl98 



AP \ ^ U 

. i — 1 



where & is the scale factor of the i* pulsar, with period 
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derivative Pi and /< is the fraction of 4ir sr of the sky cov- 
ered by the radio beam. In the discussion below, we quote 
results for 'potentially observable pulsars' (i.e. fi — 1) and 
assuming various beaming models for which fi < 1. The 
uncertainty in the current for each period bin is dominated 
by the statistical erro r in the scale factors. These ar e calcu- 
lated, again following I Vivekanand fc Naravar] Jl98lf) . as the 
square root of the sum of the squares of the scale factors in 
each period bin. 

A histogram of pulsar current versus period is shown 
for the population as a whole in Fig. [S] from which we infer 
a birth rate (i.e. maximum current value) of potentially ob- 
servable pulsars with L > 0.1 mjy kpc 2 to be 0.34±0.05 pul- 
sars per century. As in earlier analysis, we find here no 
requirement for the birth (or "injection") of pulsars with 
periods ~ 0.5 s. Following these authors, we have also ex- 
amined the birth rate contribution as a function of inferred 
dipole magnet ic field strength, B. Us ing the same group- 
ing criteria as IVranesevic et all i2004t) , we find birth rates 
of 0.04 ± 0.02, 0.14 ± 0.07 and 0.16 ± 0.04 pulsars per cen- 
tury respectively for the low (B < 9 x 10 11 G), medium 
(9 x 10 11 < B < 2.5 x 10 12 G) and high (B > 2.5 x 10 12 G) 
field populations . Thes e results agree with the findings of 
IVranesevic et alJ i2004) . 

In Table |H| we use various beaming models to apply a 
period-dependent correction to the pulsar current which ac- 
counts for those pulsars whose radiation beams do not in- 
tersect with our line of sight. We also break down these 
results for three different luminosity cut-offs: 0.1, 1 and 10 
mjy kpc 2 . For the three most recently published beaming 
models considered llLvne fc Manchester! Il988t iBiggsl [l99Ci: 
iTauris fc Manchesteilll998t) the results are consistent within 
the uncert ainties derived from the s c ale fa ctor statistics. For 
the earlier iNaravan fc Vivekanand! jl983t) model, in which 
the radio beams are elongated in the meridional direction, 
the birth rates are systematically lower. Regardless of the 
beaming model adopted, this compilation suggests that the 
majority (~ 80%) of all pulsars appear to be born with lu- 
minosities L > 1 mjy kpc 2 . It should be stressed, however, 
that the pulsar sample from which these results were ob- 
tained do not include the recent discoveries made in deep 
targeted searches of s upernova remnants and X-ray point 
sources iCamilol 120041 and references therein). As was the 
case for the luminosity function discussed in H4.5.3I further 
insights into the birth luminosities of pulsars will await a 
statistical analysis of these discoveries from the targeted 
searches. 

4.7 Predictions for other pulsar surveys 

We have shown in H4.5.2l that our models have generally pro- 
vided a good match to the yields of existing pulsar surveys. 
It is therefore natural to apply these models to make some 
predictions of the number of pulsars expected in on-going 
surveys. In Table [5] we list the number of expected pulsar 
detections for two such projects: the Perseus Arm (PA) sur- 
vey being carried out at Parkes and PALFA, a pulsar survey 
with the 7-beam Arecibo L-band Feed Arra y (ALFA) sys- 
tem dCordes et all2006t lEorimer et af]l2006ft . Note that the 
predictions below refer to normal pulsars only. 

The PA survey is extending the Galactic plane coverage 
of the PMB survey to the region -160° < I < -100° and 



|6| < 5° where there are currently 17 non-recycled pulsars. 
With the exception of a data sampling interval of 125/is, 
the PA survey uses identical observing parameters and data 
acquisition and processing schemes to those implemented in 
PMB survey. Since the PA survey is targeting pulsars out- 
side the Solar circle, the expected yield is strongly depen- 
dent on the radial distribution in the outer Galaxy. This can 
be seen by comparing the number of detections for model 
S (62) and C' (32). Upon completion, it is expected that 
this survey will provide excellent constraints on the radial 
distribution for R > 8.5 kpc. 

The PALFA survey is currently observing two regions 
of the Galactic plane visible from Arecibo with |6| < 5° at 
an observing frequency of 1.4 GHz. The inner Galaxy survey 
(PALFAi) covers the region 32° < I < 77° while the anti- 
centre survey (PALFAa) covers the region 168° < I < 214°. 
A detailed description of th e survey parameters and initial 
discoveries can be found in ICordes et al.1 ll2006h . Prior to 
these surveys, the number of normal pulsars known in these 
ranges was 200 for the inner Galaxy and 6 for the anticentre. 
Models S and C therefore predict of order 350 new discover- 
ies in the inner Ga laxy. While th i s is sig nificantly lower than 
estimates made bv lCordes et al.l (|2006) , it is consistent with 
recent independent work carried out by Faucher-Giguere & 
Kaspi (private communication). As for the PA survey, the 
lower expected PALFA yields in the anticentre region (of 
order 50 or 25 discoveries for models S and C) depend sen- 
sitively on the radial density function in the outer Galaxy. 



5 CONCLUSIONS 

We have presented discovery and follow-up parameters for 
142 pulsars discovered in the Parkes multibeam survey. Us- 
ing a sample of 1008 normal pulsars from this survey and 
the Parkes high latitude survey we have carried out Monte 
Carlo simulations to investigate various aspects of the Galac- 
tic population. Our main conclusions are as follows: 

• The derived Galactic distribution of pulsars and as- 
sumed distribution of free electrons in the Galaxy are 
strongly coupled. Bas ed on the electron density model of 
ICordes fc Laziol i2002T) , we infer a radial distribution for the 
pulsar population which peaks at 4 kpc from the Galactic 
centre. Although further surveys at higher frequencies will 
be useful, it is more important to have independent con- 
straints on the distances to currently known pulsars in this 
region to confirm the proposed radial distribution. 

• The z distribution as inferred from our simulations is 
also stron gly influenced by the el ectron density model used. 
Using the lCordes fc Laziol j2002l) model in our simulations, 
the best-fitting model population under-predicts the number 
of pulsars found in the Swinburne intermediate and high lat- 
itude surveys. To overcome this bias, we have increased the 
scale height of the model population to 330 pc. We believe 
that this apparent inconsistency will provide a very useful 
constraint for future revisions of the Galactic distribution of 
free electrons. 

• The luminosity distribution we derive follows a simple 
power law with a slope of dlog N/dlog L ~ —0.8. While this 
is somewhat flatter than previous estimates, our luminos- 
ity function is consistent with that derived independently 
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bv lFaucher-Gieuere fc Kaspi (2006). The shape of the lumi- 
nosity function below 0.1 mjy kpc 2 is not well constrained 
by our analysis. Although 6 pulsars with L < 0.1 mjy kpc 2 
are known in our current sample, their existence can be ex- 
plained by distance uncertainties and a luminosity function 
with a lower bound at 0.1 mjy kpc 2 . We stress again here 
that an analysis of recent discoveries of faint young pulsars 
in targeted searches will provide further constraints on the 
luminosity function of pulsars below 0.1 mjy kpc 2 . 

• Applying a pulsar current analysis to the observed sam- 
ple and a modified version of model C, we find the birth rate 
of the potentially observable po pulation to be 0.34±0.05 pul - 
sars per century. Applying the lTauris fc Manchester! (I1998T) 
beaming model to account for unbeamed pulsars, we find a 
birth rate of 1.4 ± 0.2 pulsars per century for luminosities 
above 0.1 mjy kpc 2 . There is no evidence for any injection 
of pulsars with birth peri ods ~ 0.5 s into the popu lation. We 
agree with the findings of IVranesevic et aD i2004T) that most 
of the birth rate comes from high magnetic field pulsars. 

• Predictions for current surveys of the outer Galaxy at 
Parkes and Arecibo depend sensitively on the form of the 
radial distribution outside the Solar circle. Although this 
distribution is poorly known at present these surveys are ex- 
pected to provide further constraints in the near future. The 
Arecibo multibeam survey of the inner Galaxy is expected to 
yield of order 350 new discoveries. We expect these surveys 
to provide further constraints on the population parameters 
derived in this work. 

The Parkes multibeam surveys provide a superb 
database with which to study various aspects of the pul- 
sar population. The analysis presented here has attempted 
to be as straightforward as possible with the minimum of 
assumptions made. F urther more detailed stud i es suc h as 
presented recently bv iFaucher-Gieuere fc Kaspil <|2006T) are 
warranted and will certainly provide further insights into 
the underlying population properties discussed here. 
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Figure 1. Mean 1400 MHz pulse profiles for 142 pulsars discovered in this phase of the multibeam survey. The highest point in the p rofile 
is pla ced at phase 0.3. Profiles have been compensated for the effects of the high-pass filter in the digitization system iManchcstcr "et* al J 
2001). For each profile, the pulsar Jname, pulse period (s) and dispersion measure (cm -3 pc) are given. The small horizontal bar under 
the period indicates the effective resolution of the profile by adding the bin size to the effects of interstellar dispersion in quadrature. 
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Figure 1. - continued 
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Figure 3. Distribution of observed beam positions in Galactic coordinates. Note that, for clarity, the size of each dot is much smaller 
than a beam width. The positions are spaced by the full-width half maximum and form an essentially uniform coverage of the region 
— 100° < I < 50°. At lower longitudes, i.e. outside the main survey region, the coverage is incomplete. A file containing the exact sky 
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Figure 4. Distribution of observed pulse widths as a function of pulse period for the observed sample (top) and a model sample (bottom). 
The model pulsars were generated assuming an empirically derived correlation given in Equation 151 
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Figure 5. Underlying (dotted lines) and derived (points) radial density functions run for two tests of the Galactic distribution deter- 
mination. The top panel shows an input population with a Gaussian radial profile, while the lower panel shows a population generated 
using Yusifov & Kiiciik's (2004) radial density profile (see section T4.3I . 
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Figure 6. Observed number distribution from our input sample (upper panels) and derived distributions for model S (lower panels) 
for the parameters: (a) p(R); (b) z; (c) L; (d) P. The solid curves are the smooth analytic functions fitted to the data (see <I4.4I . The 
dotted curves show: (a) the assumed radial density function of free electrons used ; (b) an exponential z dist r ibutio n with a scale height 
of 350 pc; (c) a log- normal fit to the opti mal pulsar population model derived bv lFaucher-Giguere &: Kaspil J2OO6); (d) a parent period 
distribution used bv lKolonko et aT] i2004T) in a study of pulse-width statistics. 
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Figure 7. Observed number distribution from our input sample (top panels) and derived distributions for model C (lower panels) for 
the parameters: (a) p{R)\ (b) z; (c) L; (d) P. Solid and dotted curves are as described in the caption to Fig. 151 
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Figure 8. Histograms of pulsar current as a function o f peri od showing the distribution for potentially observable pulsars (top) and after 
applying the beaming model of lTauris &: Manchester! Jl998t below) with luminosities greater than 0.1 mjy kpc 2 . The error bars shown 
are statistical estimates based on the sum of the squares of the individual scale factors in each bin jVivekanand & Naravar]ll98lfl . 
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Table 1. Positions, flux densities and pulse widths for 142 pulsars discovered in this phase of the Parkes multibeam pulsar survey. Radial 
angular distances are given in units of beam radii. Pulse widths at 10 percent of the peak are given only for high signal-to-noise profiles. 
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Table 1. - continued 
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6 


0.51 


9.7 


0.06 


17 


42 


1836-1324 


18 


36 


52.287(9) 


-13:24 


33.6(5) 


19.35 


-2.91 


4 


3.83 


20.0 


0.10 


2.4 


41 
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Table 1. - continued 



PSR J 


K.A. (J2000) 
(h m s) 


Dec. (J2000) 
(° ' ") 


I 

(°) 


b 

(°) 


Beam 


Radial 
Distance 


O /T\T 

o/N 


i>1400 
(mjy) 


w 50 

(ms) 


Wio 
(ms) 


1837+0053 


18:37:28.45(7) 


+00:53:13(3) 


32.13 


3.51 


7 


1.17 


18.7 


0.34 


49 


108 


1839-0436 


18:39:00.064(9) 


-04:36:57.5(7) 


27.41 


0.65 


8 


0.65 


10.9 


0.23 


6 


12 


1839-1238 


18:39:43.708(19) 


-12:38:40.7(18) 


20.35 


-3.18 


4 


3.73 


75.7 


0.37 


30 


53 


1840+0214 


18:40:34.087(19) 


+02:14:54.7(6) 


33.70 


3.44 


3 


0.31 


13.8 


0.07 


10 


22 


1840-0840 


18:40:51.9(4) 


-08:40:29(15) 


24.01 


-1.62 


2 


0.03 


61.7 


1.00 


188 


343 


1840-1207 


18:40:53.782(8) 


-12:07:32.6(6) 


20.94 


-3.20 


8 


0.77 


24.2 


0.22 


15 


22 


1841+0130 


18:41:17.761(8) 


+01:30:16.9(3) 


33.12 


2.94 


12 


1.06 


11.4 


0.06 


6 


11 


1841-1404 


18:41:34.3(4) 


-14:04:19(41) 


19.28 


-4.23 


7 


0.15 


19.0 


0.18 


17 


122 


1842+0257 


18:42:30.93(4) 


+02:57:58.8(14) 


34.56 


3.34 


4 


1.06 


55.0 


0.26 


52 


85 


1842+0358 


18:42:17.022(6) 


+03:58:35.3(3) 


35.43 


3.85 


6 


0.63 


20.2 


0.09 


4 


8 


1842+0638 


18:42:10.002(17) 


+06:38:15.0(7) 


37.80 


5.07 


10 


0.98 


16.1 


0.21 


13 


20 


1843-1507 


18:43:33.607(7) 


-15:07:03.0(8) 


18.56 


-5.13 


9 


0.63 


46.1 


0.17 


7 


18 


1845+0623 


18:45:08.59(3) 


+06:23:57.6(9) 


37.92 


4.31 


1 


0.53 


39.0 


0.33 


13 




1845-0826 


18:45:05.693(19) 


-08:26:39.7(12) 


24.69 


-2.44 


10 


0.57 


15.4 


0.33 


13 


31 


1845-1114 


18:45:45.778(6) 


-11:14:11.0(5) 


22.28 


-3.86 


2 


0.27 


41.9 


0.52 


4 


7 


1845-1351 


18:45:11.50(3) 


-13:51:54.5(20) 


19.86 


-4.92 


2 


0.61 


26.2 


0.33 


63 


84 


1846-0749 


18:46:19.744(5) 


-07:49:21.4(3) 


25.39 


-2.43 


12 


1.49 


15.9 


0.35 


7 


15 


1848-0601 


18:48:20.349(14) 


-06:01:07.7(8) 


27.22 


-2.05 


4 


0.74 


18.3 


0.24 


9 


22 


1848-1150 


18:48:11.81(3) 


-11:50:09.2(18) 


22.01 


-4.66 


10 


0.40 


34.6 


0.21 


31 


42 


1849+0409 


18:49:03.471(13) 


+04:09:42.3(5) 


36.37 


2.42 


4 


0.23 


33.1 


0.10 


9 


17 


1851-0114 


18:51:16.41(3) 


-01:14:13.4(13) 


31.81 


-0.53 


2 


0.87 


18.5 


0.28 


22 


44 


1853+0853 


18:53:22.07(11) 


+08:53:17(3) 


41.07 


3.61 


5 


0.35 


19.6 


0.11 


50 


90 


1853+1303 


18:53:57.31827(8) 


+13:03:44.0884(17) 


44.87 


5.37 


1 


0.85 


15.1 


0.4 


0.59 


2.19 


1856-0526 


18:56:21.977(15) 


-05:26:56.8(8) 


28.64 


-3.57 


13 


0.86 


35.1 


0.40 


28 


39 


1901-0315 


19:01:16.33(8) 


-03:15:14.3(20) 


31.15 


-3.67 


2 


0.96 


11.2 


0.09 


10 


20 


1901+0435 


19:01:32.2(3) 


+04:35:23(10) 


38.17 


-0.15 


5 


0.43 


43.8 


0.27 


315 




1901+0621 


19:01:05.96(9) 


+06:21:19(4) 


39.69 


0.76 


10 


0.76 


14.2 


0.47 


43 




1903-0258 


19:03:30.343(11) 


-02:58:15.6(7) 


31.66 


-4.04 


9 


0.37 


22.3 


0.14 


7 




1903+0925 


19:03:55.18(20) 


+09:25:55(6) 


42.74 


1.54 


5 


1.19 


21.7 


0.20 


89 




1904-0150 


19:04:55.555(10) 


-01:50:32.5(5) 


32.83 


-3.84 


7 


0.33 


23.0 


0.09 


5 


13 


1906+0414 


19:06:57.793(7) 


+04:14:29.44(13) 


38.48 


-1.51 


6 


0.79 


12.7 


0.23 


20 




1907+1149 


19:07:37.78(6) 


+11:49:37.1(17) 


45.29 


1.83 


11 


0.55 


25.3 


0.13 


15 


30 


1910-0112 


19:10:15.7(4) 


-01:12:06(14) 


34.01 


-4.74 


10 


0.66 


11.0 


0.07 


50 


89 


1910+1256 


19:10:09.70041(6) 


+12:56:25.5276(6) 


46.56 


1.80 


12 


0.16 


17.2 


0.5 


0.14 


0.43 


1911+1347 


19:11:55.2047(3) 


+13:47:34.411(7) 


47.52 


1.81 


1 


0.31 


36.3 


0.08 


0.15 


0.44 


1913+0904 


19:13:21.061(8) 


+09:04:45.4(4) 


43.50 


-0.68 


4 


0.54 


33.3 


0.07 


3 


6 


1914+0219 


19:14:23.794(8) 


+02:19:47.8(3) 


37.63 


-4.04 


13 


0.62 


172.0 


0.75 


16 


24 


1915+0227 


19:15:02.122(5) 


+02:27:47.77(16) 


37.82 


-4.12 


9 


0.69 


56.3 


0.55 


7 


27 


1915+1410 


19:15:35.416(12) 


+14:10:51.5(3) 


48.27 


1.20 


9 


0.47 


17.7 


0.10 


15 


24 


1921+0812 


19:21:47.704(5) 


+08:12:51.86(10) 


43.71 


-2.93 


11 


0.64 


30.0 


0.66 


2 


5 


1927+0911 


19:27:25.628(8) 


+09:11:05.7(3) 


45.22 


-3.70 


2 


0.78 


14.1 


0.16 


3 


8 


1941+1341 


19:41:04.867(18) 


+13:41:30.4(6) 


50.80 


-4.47 


9 


0.57 


17.7 


0.17 


9 


33 
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Table 2. Periods, period derivatives and dispersion measures for 142 pulsars discovered in this phase of the Parkes multibeam pulsar 
survey. We also give the MJD of the epoch used for period determination, the number of TOAs included in the timing solution, the MJD 
range covered and the RMS of the post-fit timing residuals. 



PSR J 


Period, P 

(s) 


P 

(10~ 15 ) 


Epoch 
(MJD) 


Ntoa 


Data Span 
(MJD) 


Residual 

(ms) 


DM 

(cm pc) 


0733- 


-2345 


1 


.79624970117(10) 


5.704(10) 


53014.0 


27 


52574-53454 


1 


.8 


194(3) 


0932- 


-5327 


4 


.3921587127(6) 


8.37(8) 


53099.0 


21 


52772-53426 


4. 


.0 


122(11) 


1001 


-5939 


7 


.7336402649(8) 


59.9(3) 


53069.0 


39 


52826-53312 


2. 


.0 


113(3) 


1020- 


-6026 





.140479910179(8) 


6.7386(20) 


53128.0 


31 


52854-53403 


1 


.4 


445.0(20) 


1107- 


-5907 





.252773323418(6) 


0.009(10) 


53089.0 


53 


52698-53480 





.8 


40.2(11) 


1125- 


-6014 





,0026303807397848(8) 


00000401 (9) 


53141.0 


36 


52827-53454 





,006 


52 951 


1148- 


-5725 


3 


.5599372320(7) 


11.1(13) 


52837.0 


20 


52569-53104 


2 


,6 


174(5) 


1216- 


-6410 





.003539375658423(3) 


0.00000162(19) 


53391.0 


40 


53051-53732 





,011 


47.40(3) 


1308- 


-5844 





.4646997488(7) 


8.639(5) 


51560.0 


24 


52692-53104 





.4 


205.6(11) 


1355- 


-5747 


2 


,0386738254(4) 


0.71(9) 


53513.0 


31 


53283-53742 


3 


3 


229(3) 


1357- 


-6429 





.16610832750(1) 


360.184(1) 


52921.0 


28 


52738-53104 





.9 


127.2(5) 


1405- 


-5641 





.61757468514(4) 


1.198(12) 


52898.0 


21 


52692-53104 





,6 


273(3) 


1439- 


-5501 





.028634888190455(12) 


0.0001418(10) 


53200.0 


170 


52915-53716 





.025 


14.56(5) 


1519- 


-5734 





.51875778474(3) 


4.273(5) 


53105.0 


32 


52807-53403 


1 


.8 


664(3) 


1519- 


-6106 


2 


.15430706064(14) 


8.372(19) 


53048.0 


36 


52695-53402 


1 


.8 


221(3) 


1538- 


-5732 





.341212811806(6) 


4.5498(6) 


51803.0 


40 


51300-52306 


1 


.0 


152.7(8) 


1558- 


-5756 


1 


.12234168628(3) 


186.406(4) 


53139.0 


41 


52855-53422 





.9 


127.8(10) 


1600- 


-5916 


1 


.2476661740(3) 


0.89(7) 


53116.0 


31 


52807-53425 


6 


.7 


177(4) 


1604- 


-4718 





.527465874211(20) 


4.961(9) 


52897.0 


20 


52695-53100 





.5 


52.0(16) 


1611 


-5847 





.3545503166464(18) 


0.002(4) 


53071.0 


33 


52739-53403 





2 


79.9(13) 


1617- 


-4608 





.567080221721(10) 


16.4665(10) 


53129.0 


34 


52855-53403 





.6 


321.5(12) 


1622- 


-4347 





.457681471040(9) 


5.051(4) 


52897.0 


19 


52695-53100 





.2 


177(3) 


1624- 


-4411 





.233164076390(4) 


0.7992(8) 


53058.0 


38 


52695-53422 





.6 


139.4(14) 


1624- 


-4613 





.8712426056(4) 


<0.12 


53500.0 


20 


53284-53716 


4 


2 


224.2(19) 


1627- 


-4706 





.140745806053(6) 


1.7305(6) 


53165.0 


62 


52807-53523 


1 


.7 


456.1(18) 


1631- 


-4155 





.55124125377(4) 


0.0088(73) 


53082.0 


38 


52739-53425 


1 


.8 


235(4) 


1634- 


-5107 





.50735578252(6) 


1.573(3) 


52860.0 


49 


52294-53426 


3 


.4 


372.8(20) 


1637- 


-4450 





.252869846841(20) 


0.5755(20) 


53283.0 


28 


53078-53489 





.4 


470.7(8) 


1638- 


-3815 





.69826060054(20) 


0.077(3) 


52135.0 


22 


52601-53107 





.7 


238.0(11) 


1638- 


-4725 





.7639335(3) 


4 8C91 


52651.0 


109 


51856-53445 


5. 


.3 


552.1(14) 


1643- 


-4505 





.237383081875(17) 


31.833(10) 


52923.0 


17 


52739-53107 





.9 


484(3) 


1646- 


-5123 





,530075207410(20) 


2.099(12) 


52906.0 


17 


52739-53100 


1 


.0 


279(3) 


1649- 


-3805 





.262025692969(7) 


0.0367(14) 


52885.0 


21 


52695-53100 





.3 


213.8(5) 


1650- 


-4921 





.1563993674972(10) 


1.81873(20) 


53253.0 


45 


52983-53523 





.2 


229.9(5) 


1653- 


-4030 


1 


.01937155664(16) 


44(41 


53139.0 


36 


52855-53422 


5 


.8 


425(8) 


1654- 


-3710 





.93916547876(3) 


0.735(3) 


53131.0 


36 


52859-53403 


1 


.3 


302.0(12) 


1656- 


-3621 





,730133707454(16) 


1.2752(16) 


53123.0 


39 


52824-53422 


1. 


.4 


229(3) 


1700- 


-3611 


1 


,49409061528(16) 


4.324(16) 


52828.0 


27 


52584-53100 


1 


.2 


232.7(12) 


1700- 


-4012 





.283791853785(11) 


0.099(4) 


53105.0 


35 


52808-53403 


1 


.7 


385(4) 


1711- 


-4322 





,1026182883472(6) 


0.02666(5) 


52433.0 


98 


51412-53454 





.8 


191.5(7) 


1715- 


-3859 





.92810750775(3) 


4.397(7) 


53097.0 


35 


52772-53422 


6 


.4 


817(5) 


1717- 


-3953 


1 


.08552061541(18) 


0.033(19) 


53313.0 


38 


52884-53742 


9 


.6 


466(8) 


1718- 


-4539 





,59047278772(3) 


7.507(7) 


53105.0 


29 


52808-53403 


1 


.0 


254(6) 


1722- 


-4400 





.2185540958272(15) 


0.3738(4) 


53058.0 


36 


52695-53422 





.3 


219.3(5) 


1724- 


-4500 


1 


.3091087873(3) 


3.51(12) 


53255.0 


24 


53052-53458 


3 


.3 


182(3) 


1725- 


-2852 


1. 


.257787625(10) 


1.99(15) 


52197.0 


20 


52692-53103 


3 


.7 


161(10) 


1728- 


-4028 





.866342509(4) 


0.23(6) 


52077.0 


20 


52692-53103 


2 


.4 


231(10) 


1730- 


-2900 


1 


.5384267263(3) 


8.85(8) 


53128.0 


26 


52854-53402 


2. 


.9 


289(5) 


1732- 


-3426 





,33288667881(4) 


0.402(9) 


52359.0 


23 


52146-52572 


1 


.2 


513.5(18) 


1733- 


-2533 





,65979414168(4) 


2.1(5) 


53070.0 


32 


52739-53402 


2 


.5 


242(5) 
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Table 2. - continued 



PSR J 


Period, P 


P 


Epoch 


Ntoa 


Data Span 


Residual 


DM 






(*} 


v iu ; 


(m Tr>~i 




(MJD) 


(ms) 




1733- 


-2837 


0.76818455674(6) 


1.108(8) 


53128.0 


22 


52854-53402 


1.7 


225(5) 


1734- 


-2415 


0.61252373999(12) 


1.051(12) 


53164.0 


56 


52856-53471 


6.4 


126.3(7) 


1736- 


-245 / 


2.64222343904(14) 


3.43(3) 


53095.0 


50 


52718-53471 


3.5 


170(4) 


1738- 


-2330 


1.97884743163(4) 


8.559(4) 


r" one r\ 

53206.0 


33 


52695-53716 


1.0 


99.3(18) 


1739- 


-3951 


0.341772476799(10) 


0.02(5) 


52897.0 


17 


52692-53103 


0.4 


24.6(4) 


1741- 


-2719 


0.346796929142(19) 


0.081(4) 


rr n n 

52836.0 


23 


52569-53103 


1.2 


361.9(4) 


1742- 


-3957 


1.0163491951(4) 


0.082(29) 


5olz4.U 


29 


coonc /loo 

d2#zd-qo422 


7.9 


186(8) 


1 T A A 

1 / 44- 


-3922 


0.172444360995(2) 


0.00155(12) 


52530.0 


36 


51953-53107 


0.3 


148.1(7) 


1750- 


-2043 


5.639047079(3) 


7.9(16) 


52923.0 


18 


52739-53107 


6.2 


oon / *7\ 

239(7) 


1751- 


-2857 


0.0039148731963690(6) 


0.00001126(4) 


52560.0 


168 


51808-53312 


0.029 


42.808(20) 


1753- 


-1914 


0.0629548889725(12) 


0.00202(12 ) 


53150.0 


62 


52828-53471 


0.8 


105.3(3) 


1754- 


-3510 


0.392703893399(3) 


0.7816(5) 


53097.0 


43 


52733-53461 


0.4 


82.3(3) 


1755- 


-1650 


0.73374443871(3) 


0.686(3) 


53165.0 


55 


52860-53471 


1.1 


159.9(10) 


1756- 


-2251 


0.02846158845494(2) 


0.0010171(2) 


52086.0 


1597 


50997-53176 


0.042 


121.18(2) 


1758- 


-2846 


0.76670638621(3) 


0.094(18) 


53253.0 


39 


53063-53443 


1.3 


66.6(3) 


i / 59- 


-1736 


0.79845156695(8) 


0.255(20) 


53156.0 


18 


52854-53458 


1.7 


206(3) 


1801- 


-1417 


0.0036250966601209(18) 


0.00000531(11) 


53156.0 


89 


52569-53743 


0.021 


57.21(4) 


1801- 


-3458 


1.3856036231(3) 


0.4(4) 


53112.0 


27 


52827-53397 


4.5 


146(6) 


1802- 


-2124 


0.012647593582763(5) 


0.000072(1) 


52855.0 


30 


52605-53105 


0.007 


149.6(1) 


1805- 


-1504 


1.1812692298(3) 


0.27(3) 


53035.0 


58 


52599-53471 


14.2 


225(3) 


1808- 


-1020 


0.596993202504(4) 


0.7725(4) 


53035.0 


63 


52599-53471 


0.5 


225.3(8) 


1808- 


1726 


0.24103454795(9) 


<0.012 


53108.0 


67 


52744-53471 


23.2 


536(7) 


1808- 


-2701 


2.45788177758(12) 


65.818(11) 


53112.0 


56 


52751-53472 


2.6 


95(4) 


1811- 


-2439 


0.415812941402(11) 


0.2972(11) 


53165.0 


54 


52859-53471 


0.9 


172.0(5) 


1812- 


-2526 


0.315835048197(8) 


0.1775(9) 


53171.0 


54 


52871-53471 


1.1 


361.4(4) 


1814- 


-0618 


1.37786849563(14) 


0.292(16) 


53095.0 


66 


52718-53471 


5.7 


168(6) 


1816- 


-0755 


0.2176426911444(15) 


6.48032(15) 


53035.0 


60 


52599-53471 


0.5 


116.8(4) 


1817- 


-0743 


0.438095346909(18) 


<0.0083 


52887.0 


13 


52667-53107 


0.7 
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Table 2. - continued 
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Table 3. Derived parameters for 142 pulsars discovered in this phase of the Parkes multibeam pulsar survey. Listed are the base-10 
logarithms of characteristic age (yr), the surface dipole magnetic field strength (G), the loss in rotational energy (erg s _1 ), the DM-derived 
distance using the Taylor & Cordes (1993) model, -D^c (kpc), and the Cordes & Lazio (2001) model, Dql (kpc), the corresponding z 
heights, ztc an d Z CL (both in kpc), and the inferred radio luminosity at 1400 MHz for each distance estimate L^c an d Lql (both in 
mjy kpc 2 ). 
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Table 3. - continued 
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0.16 


34.8 


24.0 


1808-2701 


5.77 


13.11 


32.23 


2.43 


1.71 


-0.14 


-0.10 


0.9 


0.4 


1811-2439 


7.35 


11.55 


32.20 


3.79 


3.62 


-0.20 


-0.19 


3.7 


3.4 


1812-2526 


7.45 


11.38 


32.34 


13.76 


8.00 


-0.83 


-0.48 


34.1 


11.5 


1814-0618 


7.87 


11.81 


30.64 


5.86 


4.22 


0.53 


0.38 


19.9 


10.3 


1816-0755 


5.73 


12.08 


34.40 


3.24 


2.78 


0.23 


0.20 


1.8 


1.3 


1817-0743 


> 8.92 


< 10.79 


< 30.59 


0.77 


0.80 


0.05 


0.05 


0.1 


0.2 


1819-1318 


7.60 


11.98 


30.83 


1.52 


1.15 


0.02 


0.02 


0.4 


0.2 


1820-0509 


6.76 


11.75 


32.98 


2.84 


2.51 


0.22 


0.20 


2.0 


1.6 


1821-0256 


8.25 


11.10 


31.32 


2.55 


2.36 


0.24 


0.22 


1.2 


1.1 


1822-0848 


7.99 


11.53 


30.96 


4.24 


3.77 


0.17 


0.15 


0.7 


0.6 


1824-0127 


7.00 


12.50 


31.00 


2.21 


1.88 


0.20 


0.17 


2.9 


2.1 


1824-2233 


7.79 


11.78 


30.88 


4.22 


3.78 


-0.33 


-0.29 


3.9 


3.1 


1824-2328 


7.14 


12.21 


31.30 


5.61 


4.55 


-0.48 


-0.39 


10.1 


6.6 


1827-0750 


6.44 


11.82 


33.49 


7.45 


6.04 


0.23 


0.19 


77.7 


51.1 


1828-2119 


6.81 


11.91 


32.57 


12.06 


6.77 


-0.99 


-0.56 


55.3 


17.4 


1829+0000 


6.78 


11.51 


33.41 


3.40 


3.22 


0.29 


0.27 


5.0 


4.5 


1830-0052 


7.36 


11.46 


32.36 


7.37 


5.68 


0.55 


0.42 


2.2 


1.3 


1830—0131 


6.06 


11.76 


34.36 


2.69 


2.66 


0.19 


0.19 


2.5 


2.5 


1830-1414 


8.20 


11.39 


30.81 


7.80 


6.26 


-0.25 


-0.20 


6.1 


3.9 


1831-0952 


5.11 


11.88 


36.04 


4.32 


4.04 


-0.01 


-0.01 


6.2 


5.4 


1832+0029 


6.75 


11.96 


32.59 


1.45 


1.32 


0.11 


0.10 


0.3 


0.2 


1834-0031 


7.06 


11.59 


32.69 


4.47 


3.96 


0.27 


0.24 


3.4 


2.7 


1834-1855 


7.08 


12.23 


31.38 


5.95 


4.62 


-0.52 


-0.40 


17.0 


10.2 


1835-0349 


6.64 


12.21 


32.30 


5.27 


5.59 


0.17 


0.18 


4.4 


5.0 


1835-0944 


5.72 


11.91 


34.75 


4.63 


4.37 


-0.08 


-0.08 


8.8 


7.8 


1835-0946 


8.14 


11.11 


31.49 


4.02 


3.61 


-0.07 


-0.06 


2.9 


2.3 


1835-1548 


6.79 


12.04 


32.36 


13.52 


7.51 


-0.89 


-0.50 


11.0 


3.4 


1836-1324 


6.44 


11.64 


33.86 


3.87 


3.48 


-0.20 


-0.18 


1.5 


1.2 
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Table 3. - continued 



PSR J log[r c ] log[B] log[£] D T c D C l ztc zcl L T c L C l 
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71 
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7. 


01 
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-0 
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-0 


42 
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1846-0749 
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11. 
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33.08 


9 
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7. 
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42 
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33 
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17.9 
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7. 
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.14 


32.00 


7. 
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Table 4. Summary of the 13 binary pulsars discovered in this phase of the Parkes multibeam survey. For each system, we give the spin 
period (P), orbital period (Pb), projected semi-major axis (x) and orbital eccentricity (e). The right-hand column gives the source of the 
more detailed parameters which either refer to external papers or additional tables below. 



PSR 


P 


P b 

( H avs ) 


X 

(It sec) 


e 








Further details 


J1125-6014 


2.63 


8.75 


8.34 


7. 


9 x 


10" 


-7 


Table El 


J1216-6410 


3.54 


4.04 


2.94 


6. 


.8 x 


io- 


-6 


Table El 


J1439-5501 


28.6 


2.12 


9.84 


5. 


.0 x 


io- 


-5 


Table 


J1638-4725 


764 


1940 


2380 


0. 


95 






Lyne et al. (in preparation) 


J1711-4322 


103 


922 


140 


2. 


.4 x 


io- 


-3 


Table El 


J1744-3922 


172 


0.191 


0.212 


1. 


.3 x 


io- 


-3 


Faulkner et al. f 20041 


J1751-2857 


3.92 


111 


32.5 


1. 


.3 x 


io- 


-4 


Stairs et al. (20051 


J1756-2251 


28.5 


0.320 


2.76 


0. 


18 






Faulkner et al. (2005) 


J1802-2124 


12.6 


0.699 


3.72 


3. 


.2 x 


io- 


-6 


Faulkner et al. ("20041 


J1822-0848 


835 


287 


97.8 


0. 


059 






Table El 


J1841+0130 


29.8 


10.5 


3.50 


8. 


.2 x 


io- 


-5 


Table El 


J1853+1303 


4.09 


116 


40.8 


2. 


.4 x 


io- 


-5 


Stairs et al. ( 20051 


J1910+1256 


4.98 


58.5 


21.1 


2. 


.3 x 


10" 


•4 


Stairs et al. ("20051 



Table 5. Orbital parameters for the two binary pulsars obtained using the Blandford & Teukolsky (1976) binary model. The minimum 
companion mass is calculated by assuming an inclination angle of 90° and a neutron star mass of 1.35 Mq. Figures in parentheses 
represent 1-<t uncertainties in the least significant digit (s). 



PSR 


J1711-4322 


J1822-0848 


Orbital period (d) 


922.4707(7) 


286.8303(14) 


Projected semi-major axis of orbit (It sec) 


139.6245(6) 


97.7279(20) 


Eccentricity 


0.002375(6) 


0.058962(9) 


Epoch of periastron (MJD) 


50208.8(3) 


52672.848(9) 


Longitude of periastron (degrees) 


293.75(12) 


272.583(10) 


Mass function (Mq) 


0.00343446(2) 


0.0121811(3) 


Minimum companion mass (Mq) 


0.20 


0.32 



Table 6. Orbital parameters for the four binary pulsars obtained using the 'ELL1' binary timing model iLange et al.ll200l|l . where the 
first and second Laplace-Lagrange parameters are defined by e cos w and e sin u> respectively (for an orbital eccentricity e and longitude of 
periastron ui). The minimum companion mass is calculated by assuming an inclination angle of 90° and a neutron star mass of 1.35 Mq. 
Figures in parentheses represent 1-cr uncertainties in the least significant digit (s). 



PSR 


J1125-6014 


J1216-6410 


J1439-5501 


J1841+0130 


Orbital period (d) 


8.75260353(5) 


4.03672718(6) 


2.117942520(3) 


10.471626(5) 


Projected semi-major axis of orbit (It sec) 


8.339198(5) 


2.937091(7) 


9.833404(7) 


3.50409(18) 


First Laplace-Lagrange parameter 


-0.0000008(13) 


0.0000004(60) 


-0.0000495(15) 


-0.0000033(10) 


Second Laplace-Lagrange parameter 


0.00000005(120) 


-0.000007(6) 


0.0000059(13) 


-0.00008(10) 


Epoch of ascending node (MJD) 


53171.5856408(11) 


53055.3611070(20) 


53058.0499527(3) 


52747.55445(11) 


Mass function (Mq) 


0.008127952(7) 


0.001669461(6) 


0.2275967(2) 


0.00042129(3) 


Minimum companion mass (Mq) 


0.28 


0.16 


1.11 


0.096 
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Table 7. Summary of the analytic function fits to models S and C. The various parameters used are described in Equations UOlTHjl 
Figures in parentheses give the 1-cr uncertainty in the least significant digit. 



Parameter 


Model S fit 


Model C fit 


Unit 


A 


44(7) 


41(5) 


kpc -2 


13 


0.2(2) 


1.9(3) 




C 


1.4(6) 


5.0(6) 




D 


0.39(2) 


0.75(3) 




E 


0.33(3) 


0.18(1) 


kpc 


F 


-0.59(5) 


-0.77(7) 




G 


3.45(8) 


3.5(1) 




H 


0.51(1) 


0.52(1) 




1 


2.71(1) 


2.70(1) 




J 


-0.34(1) 


-0.34(1) 





Table 8. A comparison of the birth rates derived from the pulsar current analysis as a function of luminosity cut-off and beaming model. 
Figures in parentheses represent statistical 1-cr uncertainties in the least significant digit(s) calculated as the square root of the sum of 
the squares of the scale factors in the maximum pulsar current bin (Vivckanand & Naravan 1981). 



Beaming model Minimum luminosity (mjy kpc 2 ) 



0.1 1 10 



No correction (i.e. / = 1) 0.34(5) 0.28(4) 0.12(3) 

Tauris fc Manchester (19981 1.38(21) 1.14(18) 0.42(3) 

Biggs (19901 1.14(18) 0.94(15) 0.35(8) 

Lvne fc Manchester (19881 1.32(21) 1.09(17) 0.41(9) 

Naravan fc Vivekanand (19831 0.51(8) 0.43(7) 0.16(4) 



Table 9. A comparison of the number of detected pulsars for the real sample with several models developed in this paper for the various 
pulsar surveys. T he completed pulsar surveys are the Parkes multi beam survey (this p aper, PMB), the Swinburne intermediate latitude 
survey lEdwards et all200ll SIL), the Parkes high latitude survey iBurgav et all2006l PH), the Swinburne high latitude survey Jjacobvl 
2005, SHL). Pulsar surveys wh ich are currently on -going are the Perseus arm (PA) survey with the Parkes multibeam system and the 
PALFA surveys with Arecibo JCordes et al.ll2006l) . PALFAi refers to the inner Galaxy survey, while PALFAa refers to the anticentre 
survey. For the completed pulsar surveys, the column labeled 'Real' lists the number of pulsars detected (i.e. discoveries and redetections). 
For the on-going surveys the numbers marked with a * refer to currently known pulsars in the survey regions with flux densities higher 
than the nominal survey threshold. The remaining columns are the average number of predicted detections using each model. 



Survey Real Model S Model C Model C' 



Completed pulsar surveys 



PMB 


986 


967 


991 


984 


SIL 


156 


146 


88 


160 


PH 


32 


33 


23 


30 


SHL 


51 


52 


17 


38 




On-. 


going pulsar 


surveys 




PA 


17* 


62 


35 


32 


PALFAi 


200* 


480 


550 


530 


PALFAa 


6* 


54 


30 


30 



